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Abstract 



We investigate effective Yukawa couplings of mesons to the elementary 
fermions in noncompact lattice QED. Phe couplings are extracted from 
suitable fermion-antifermion-meson three-point functions calculated by 
Monte Carlo simulations with dynamical staggered fermions. Phe sca- 
ling behaviour is compatible with expectations from perturbation theory, 
thus indicating triviality of QED. Phe lines of constant Yukawa coupling 
are compared to flows of other quantities. Consistency is seen, at most, 
for weak coupling. 



In strongly coupled lattice QED chiral symmetry is known to be spontaneously 
broken, whereas an unbroken chiral symmetry is expected in the weak coupling 
limit [1]. This finding immediately raises questions like: Can a nontrivial continuum 
limit be defined at the corresponding critical point? If so, what is the nature of the 
associated continuum theory? Several groups have investigated these problems by 
means of extensive Monte Carlo simulations of lattice QED using various techniques 
(see, e.g., [2-9]). Yet the physical interpretation of the data is still controversial (see 
also [10]). 

Correlation functions of composite fermion-antifermion operators show evidence 
for at least two fermion-antifermion (bound) states ("meson" states): a pseudos- 
calar state denoted by P (the Goldstone boson associated with the spontaneous 
symmetry breaking) and a scalar state (S) [3, 11, 6]. In this letter we study fermion- 
antifermion-meson three-point functions in order to determine (effective) Yukawa 
couplings of P and S to the fermions. The investigation should provide additional 
information about the lattice model. In particular, if it has a trivial continuum 
limit, the couplings should tend to zero as one approaches the critical point. 

We have performed simulations with dynamical staggered fermions using the 
noncompact formulation of the gauge field action. So the total action is given by 
S = So + Sf with 

Sa = tfT, (!) 
S f = E {§ E Vvix) (x(x)e iAl4x) x(x + A) - X(x + fi)e- iA ^x(x)) + mx(x) X (x)} 

X jJ, 

(2) 

where f3 = 1/e 2 (e = bare charge) and 

^(x) = (-ir + -+^- 1 , (3) 

F^{x) = A^x) + A v {x + fi)- A^x + v)- A v (x). (4) 

The data in this paper were obtained on an L 3 S X L t lattice with L s = L t = 12. The 
boundary conditions are periodic for the gauge field A^. For the fermion field XiX 
we have chosen periodic (antiperiodic) spatial (temporal) boundary conditions. For 
more details on the simulation see ref. [6]. 

In order to calculate fermionic observables like the fermion propagator or the 
fermion-antifermion-meson three-point functions, which are the subject of this paper, 
we have to fix the gauge. We choose to work with the Landau gauge, 

Y J (A,(x)-A,(x-f 1 )) = 0, (5) 

A 1 
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which has the advantage that it can be implemented exactly in the noncompact 
formulation. Unfortunately, the zero-momentum mode of the gauge field does not 
average to zero in our ensembles. To cope with this problem we divide our data into 
blocks of, say, 20 successive configurations, on which this mode is approximately con- 
stant, and perform our analysis in each block separately taking the zero-momentum 
mode into account explicitly (cf. [6]). The averages and errors are then obtained 
from these block results. 



As we have argued in ref . [6] , the renormalized charge vanishes as one approaches 
the critical point f3 c ~ 0.19 , m = 0. Moreover, it turns out that the chiral phase tran- 
sition is reasonably well described by logarithmically improved mean field theory, 
i.e., the equation of state can be derived from an effective linear cr-model [6, 12]. Wi- 
thin such a model, the effective quartic coupling A is driven to zero in the continuum 
limit according to a (one-loop) renormalization group equation of the form 



dA 
dt 



c\ 2 ,t = In a, (6) 



where a = (xx) i s the chiral condensate and c > 0. From the solution of this 
equation one sees that A ~ |lncr| _1 as a — > 0. On the other hand, for a Yukawa 
coupling g, which we are investigating here, one finds in leading order [13] 

, c'>0, (7) 

so that 

— = -2c't + const. (8) 
9 

Hence the triviality scenario leads to the expectation that g vanishes like 

g ~ |ln(j|- 1/2 . (9) 



We shall extract our effective Yukawa couplings from three-point functions of 
the form 

G£\p, q; t , h) = J2 e 2i ^'- 2i ^(x(2f , h)M a (p - £ 0) X (2f, t )) c . (10) 

x,x r 

Here (• • -) c denotes the fermion-line connected part of the expectation value, the 
lattice vectors x,x' label spatial cubes of size 2 3 , and M a (p 0} t) is the standard 
pseudoscalar (a = P) or scalar (a = S) meson operator: 

M a (po,t) = ^ M ^{y)xiMx{M (11) 

y 
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with (fip(y) = ( — l) yi+y2+ys } (fis(y) = 1- We have calculated the three-point function 
(10) for the following values of the spatial fermion momenta p } q: 

, _> -.x ,-> -> ^ /27r_, -> . ,2tt 2tt . 

(p, q ) = (0,0), ( — e 3 , ), (— e 3 , — e 3 ) . (12) 

The spatial meson momentum is then given by p = p — q. 

As is well known, the operators (11) couple not only to the state one is interested 
in but also to its "parity partner" (and higher excitations). Having computed the 
three-point function (10) for arbitrary times t 0} t\ we can, however, project onto a 
single parity by Fourier transforming the t 0} ti dependences: 

E e_iP4tl+i,4toG! 5 a) (p>?;*o,*i) 

= \ E e-^ tl+i ^ e 2i9 "'- 2i? %(2f, h) e^-^M^p - q, t) X (2x', t )) c . 

* ti,to x,x' t 

(13) 

Since the Goldstone boson is represented by the component of Mp that "oscillates" 
in time (see, e.g., [14]), we must choose p 4 — q 4 = tt in order to study the Yukawa 
coupling of the pseudoscalar. In the case of the scalar, on the other hand, we have 
to take P4 — q 4 = 0. 



We define effective Yukawa couplings g a by comparing the Monte Carlo data 
for the three-point functions with tree-level formulas derived from an effective lat- 
tice action which describes the interaction between staggered fermions XiX an d a 
(pseudo-) scalar field $ via coupling terms of the form 

-gsJ2®( x )x( x )x( x ) (scalar), (14) 

X 

—gp^2^(x)x(x)x( x )( — ^T 1+ "' +Xi (pseudoscalar). (15) 

X 

In the tree-level formulas for the three-point functions we replace the free propa- 
gators by the full propagators calculated in the simulations. Furthermore, we in- 
clude the appropriate wave function renormalizations Zp (for the fermion field) and 
Z a (for the mesons). In this way we arrive at the following equation defining the 
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where 



momentum-dependent Yukawa couplings g a (p } q): 

E e_iP4tl+i,4to G ! 5 a) (p>?;*o,*i) 

_1 _1 _1 

= 9a(p,q)Z F (p) 2 Z F (q) 2Z a (p-q) 2 

X E j {Pi - 9i)t js(M a (q-p, t)M a (p - q, 0)) c (lg) 

a; to 27, of' 

£ ^) V^+^ £ e^*"* <X(25 + *i)x(2^, 0))}* , 

x,x' 

S = -i, e P = 1, 

(17) 

# 5 R = (-ir + ^ 3 ,M^) = i, 

and o3 runs over the 8 three- vectors with components or 1, i.e. uj labels the corners 
of a spatial unit cube. 

The wave function renormalizations Zp,Z a are obtained in the usual way by 
comparing the one-particle contributions in the full propagators with the correspon- 
ding free propagators. For example, writing the one-particle contribution to the 
meson propagator 

^(M a (-p,t)M a (p,0)) c (18) 

in the form 
we have 

Z a (p) = 2\A a (p)\ (l- e - LtE °W) sinhE a (p) (20) 

(see, e.g., ref. [15]). In the case of the fermion wave function renormalization we 
compare with the propagator in a constant background potential [6]. 

Let us now discuss the results. One first observes that g a (p, q) is essentially real: 
The imaginary part is relatively small and noisy, for special values of the momenta 
it is even exactly zero. Hence we shall disregard it and consider real quantity. 

For obtaining a final unique answer it would be desirable to send all momenta to 
zero. However, the antiperiodic temporal boundary conditions for the fermions 
would allow us only to give the mesons momentum zero by choosing p = q,p± = q± 
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for the scalar and p 4 = g 4 + tt (mod 2tt) for the pseudoscalar. Admitting for (p } q ) 
the values (12) we give |g 4 | the smallest possible value and average over q 4 = ±7r / L t : 

9s(p,q) = \ J2 9s((p,q4),(q,q4)), (21) 

g 4 =±7r/L t 

gp(p,q) = l J2 9p((p,q4 + *),(q,q4))- (22) 

g 4 =±7r/L t 

For the error we take the largest error of a single value. 

In fig. 1 we plot gp(p } q)~ 2 versus |ln<r| as suggested by eq. (8). (Note that in 
the preliminary presentation of our results in [16] we have erroneously shown the 
average of <?p((0, 54), (0, ^4))-) Although these numbers were obtained at different 
bare couplings (3 (0.17 < (3 < 0.22) and different bare masses m (0.02 < m < 0.09), 
they fall into a narrow band, which with some optimism could even be considered 
as a single straight line in accordance with the triviality scenario (cf. eq. (8)). We 
consider these data as a qualitative confirmation of the triviality hypothesis. The 
results for the coupling of the scalar shown in fig. 2 point into the same direction. 

However, several caveats have to be kept in mind. We have neglected the fermion- 
line disconnected contributions to the three-point functions and to the meson propa- 
gators, since they are difficult to calculate. These contributions could be especially 
important in the case of the scalar, which is a flavour singlet. Furthermore it is 
not clear, how our Yukawa couplings are affected by the fact that for some of our 
data points the mesons are heavier than two fermions and hence might correspond 
to resonances rather than stable particles in the continuum limit. Another pro- 
blem concerns the definition of the wave function renormalizations. Since they are 
extracted from fits of free propagators to the measured two-point functions, they 
correspond to on-shell renormalization, whereas the effective Yukawa couplings are 
defined at fixed momenta ("fixed" in lattice units, the physical scale varies from si- 
mulation point to simulation point). The weak momentum dependence of the results 
may however indicate that this problem is not too severe. 

With the effective Yukawa couplings at hand we can draw new flow diagrams, 
i.e. we can construct the lines of constant gs, gp in the (3-m plane. These are to 
be compared with the lines of constant renormalized charge ep or constant mass 
ratios [4, 6]. In fig. 3 we show the curves of constant gp(0,0) together with the 
lines of constant ep. We recall that the curves of constant ep > end on the line 
of first-order chiral phase transitions m = 0,f3 < (3 C} while only the curve ep = 
flows into the critical point [6], in accordance with triviality. The lines of constant 
<7p(0,0) look rather similar, although truly parallel flows seem to be possible only 
in the lower right hand corner, i.e. for small values of ep and gp. This suggests that 
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also in the case of the Yukawa couplings only the curve of vanishing coupling flows 
into the critical point as perturbation theory would predict. 

In fig. 4 we compare with the lines of constant nip/ nip where nip is the renor- 
malized fermion mass and nip the pseudoscalar mass. Here we observe a completely 
different flow pattern in the left part of the diagram, whereas there is a tendency 
towards parallel flows in the region of small ep. The corresponding plots for gs look 
qualitatively similar. 

In conclusion, these results are in agreement with our previous findings [4, 6] 
indicating that lines of (approximately) constant physics can only be expected for 
small couplings. To achieve renormalizability for larger couplings it would seem to 
be necessary to include further relevant operators in the action, such as four-fermion 
terms. 
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Figure 1: g P (p,q)- 2 versus ln(l/<r) for (j3,(f) = (0, ), (ge 3 , ge 3 ), (ge 3 , ) (top to 
bottom). 
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Figure 2: gs{p,q)~ 2 versus ln(l/<r) for (p, <f ) = (0, ), (ge 3 , ge 3 ), (ge 3 , ) (top to 
bottom). 
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Figure 3: Curves of constant gp(0, ) (full lines) and curves of constant renormalized 
charge (dashed lines) with gp(0, ) ranging from 2.0 (lower right-hand corner) to 3.6 
(upper left-hand corner) in steps of 0.2 and e 2 R ranging from 2.8 (lower right-hand 
corner) to 4.6 (upper left-hand corner) in steps of 0.2. 
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Figure 4: Curves of constant </p(0,0) (full lines) and curves of constant mp/mp 
(dashed lines) with gp(0, ) ranging from 2.0 (lower right-hand corner) to 3.6 (upper 
left-hand corner) in steps of 0.2 and mp/mp ranging from 0.4 (lower right-hand 
corner) to 1.6 (lower left-hand corner) in steps of 0.1. 
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